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During  the  past  three  years  we  have  investigated'  a rather  broad 
array  of  topics  that  include  microwave  and  quasi-optical  Integrated 
circuits,  coherent  optics  and  the  application  of  numerical  and  asymptotic 
techniques  to  radar  scattering  and  antenna  problems.  In  the  following 
sections  we  briefly  describe  some  of  the  important  accomplishments  in 
these  areas  and  list  the  scientific  publications  describing  the  research 
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II.  TECHNICAL  DESCRIPTION 


I.  Integrated  Circuits 

Work  on  components  for  microwave  integrated  circuits  was  continued 
from  the  past  grant  period  with  emphasis  being  placed  on  the  development 
of  efficient  techniques  for  analyzing  distributed  circuit  components.  To 
illustrate  the  application  of  the  transform  method  of  solution  which  was 
previously  developed  [1],  the  problem  of  computing  the  characteristics  of 
microstrip  resonators  was  investigated  and  the  numerical  results  were 
reported  [2]  for  a resonator  of  rectangular  shape. 

New  directions  of  research  on  integrated  circuits  were  initiated 
during  this  grant  period  with  a view  to  bridging  the  gap  between  the 
microwave  and  optical  integrated  circuits.  This  effort  led  to  the  develop- 
ment of  new  concepts  in  quasi-optical  and  millimeter  waveguides  as  well 
as  components  [3,4]  suitable  for  integrated  circuit  systems  in  the  20-90 
GHz  (or  even  wider)  frequency  range.  Various  types  of  planar  waveguide 
structures  with  different  types  of  materials  were  investigated.  The 
designs  that  showed  the  most  promise  appeared  to  be  the  "Inverted  Strip",  or 
the  IS  guidejWhich  consists  of  a dielectric  (6^)  slab  mounted  on  a 
dielectric  (G2>  guide  of  rectangular  cross  section  resting  on  a conducting 
plane.  An  extensive  study  of  the  IS  guide  was  carried  out  for  the  choice 
of  quartz  for  6^  and  teflon  for  62>and  the  results  were  found  to  be  very 
encouraging.  Near-field  measurements  were  made  to  verify  that  the  inverted 
strip  design  indeed  confines  the  fields  in  a region  away  from  the  ground 
plane,  a key  requirement  for  reducing  conductor  losses.  Techniques  for 
the  theoretical  analysis  of  open  dielectric  waveguides  were  developed  using 
the  concept  of  'equivalent  dielectric  constant'  that  allows  efficient  and 
accurate  analysis  of  the  dominant  modes  propagating  in  such  waveguides. 
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The  theoretical  results  were  correlated  with  experimental  measurements  and 
good  agreement  was  demonstrated  between  the  two. 

In  addition  to  the  waveguides  a number  of  integrated  circuit  components, 
e.g.,  directional  couplers  and  resonators, were  also  investigated,  both 
theoretically  and  experimentally  and  the  results  were  reported  [5,6]. 

During  the  next  grant  period  we  plan  to  investigate  yet  another 
form  of  the  IS  guide  [7]  which  employs  a single  dielectric  (Stycast  Hi-K) 
design,  thus  making  it  convenient  for  fabrication  using, for  example,  injection 
molding  or  etching  techniques.  Refined  methods  for  analyzing  these  wave- 
guides will  also  be  developed,  particularly  with  a view  to  accurately 
deriving  the  propagation  characteristics  in  open  dielectric  waveguides. 

In  addition  to  millimeter  waveguides,  the  problem  of  excitation 
of  optical  fibers  was  also  investigated  and  the  optimum  dimensions 
of  the  system  were  computed  [8]  of  higher-order  modes  which  are  usually 
present  in  such  waveguides. 

2.  Numerical  and  Asymptotic  Techniques  in  Electromagnetics 

Extensive  research  into  important  aspects  of  numerical  and  asymp- 
totic techniques  in  Electromagnetics  has  been  carried  out  during  the 
past  three  years  under  the  present  grant.  The  effect  of  different  testing 
functions  in  the  moment  method  solution  of  thin-wire  antenna  problems  has 
been  investigated  [9], and  some  definitive  conclusions  with  regard  to  the 
moment  method  approach  to  numerically  solving  electromagnetic  problems 
have  been  derived.  We  have  also  investigated  other  anomalous  characteristics 
of  the  moment  method,  viz.,  the  unstable  nature  of  the  matrix  when  the 
frequency  of  the  incident  wave  is  in  the  neighborhood  of  the  internal 
resonant  frequencies  of  a scatter.  Methods  for  alleviationg  this  type 
of  difficulty  were  developed  and  described  in  a journal  paper  [10]. 
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A new  method  for  solving  the  problem  of  scattering  from  electrically 
large  structures  with  application  to  radar  scattering  problems  was  developed. 
The  novel  technique, called  the  "Spectral  Theory  of  Diffraction"  [11],  not 
only  provides  a new  explanation  of  the  failure  of  the  ray  theories  at  the 
so-called  "trouble-regions,"  e.g.  shadow  boundaries  and  caustics, but  shows 
how  these  difficulties  can  be  circumvented.  The  spectral  approach  is  also 
useful  for  systematic  improvement  of  high  frequency  asymptotic  solutions 
as  has  been  demonstrated  with  a number  of  practical  examples  [11]  - [13]. 

This  feature,  which  has  not  been  available  hitherto  in  the  ray  techniques, 
opens  up  new  avenues  for  solving  complex  problems  for  which  the  methods 
based  on  the  Geometrical  Theory  of  Diffraction  are  not  sufficiently 
accurate.  Further  developments  of  the  spectral  technique  are  planned  in 
the  future. 

Another  important  problem  in  the  area  of  numerical  electromagnetics 
that  we  have  investigated  involves  a monopole  antenna  situated  on  a 
finite  ground  plane  [14].  The  current  distribution  on  the  ground  plane  and 
the  antenna  itself  have  been  calculated  both  for  the  antenna  operating  in 
the  receiving  and  transmitting  mode.  The  numerical  approach  followed  in 
this  work  differs  from  conventional  techniques  and  was  chosen  because  it 
gives  stable  solutions  for  the  current  distribution  on  structures  with 
edges,  whereas  the  singularity  of  one  of  the  current  components  introduces 
difficulties  in  the  commonly  used  intregro-dif ferential  equations. 

3.  Digital  and  Optical  Processing 

Primary  emphasis  on  research  in  this  area  was  devoted  toward  the 
problem  of  digital  imaging  of  incoherent  sources  and  to  the  problem  of 
decoding,  deconvolut ing  or  reconstructing  these  images  recorded  using  a 
coding  aperture  [15,16]. 
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Relative  advantages  of  digital  processing  techniques  over  their 


optical  counterparts  have  been  pointed  out,  particularly  from 
of  view  of  reducing  the  distortion  in  the  image. 

The  problem  of  correcting  the  image  aberration  using  a 
optical  data  processing  technique  and  a holographic  technique 
deblurring  the  image  has  been  reported. 


the  point 

coherent 
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Spectral-Domain  Approach  for  Calculating  the 
Dispersion  Characteristics  of  Microstrip  Lines 

T.VTsrO  ITOII  and  l< \J  Mi  l TUA 

Abstract — The  boundary  value  problem  associated  with  the  open 
microstrip  line  structure  is  (orrnulated  in  terms  ol  a riyorous,  hybird- 
mode  representation.  The  resulting  equations  are  subsequently 
transformed,  via  the  application  of  Galerkin's  method  in  the  spectral 
domain,  to  yield  a characteristic  equation  for  the  dispersion  proper- 
ties of  the  open  microstrip  line. 
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Excitation  of  an  optical  fiber  by  a Gaussian  beam 


M.  Mostafavi,  T.  Koh,  and  R.  Mittra 


In  thii  paper  we  present  the  results  of  an  in-depth  analysis  of  the  problem  of  Gaussian  I team  excitation  of 
optical  fibers  at  normal  incidence.  The  method  entails  the  application  of  the  continuity  conditions  at  the 
interface  located  at  the  fiber  end  using  Gaussian  beams  with  different  spot  sizes  as  our  excitation  source. 
Next,  the  simultaneous  set  of  equations  obtained  is  numerically  solved  for  the  reflected,  guided,  and  ra- 
diated modal  coefficients  and  their  corresponding  powers. 


The  Effect  of  Different  Testing  Functions  in  the  Moment 
Method  Solution  of  Thin-Wire  Antenna  Problems 

C.  A.  KLEIN  and  It.  M1TTHA 

Abstract — The  use  of  piecewise  sinusoids  for  expansion  func- 
tions »nd  rectangular  pulses  for  testing  functions  is  described  in 
the  application  of  the  method  of  moments  to  thin-wire  antennas 
and  scattcrers.  This  choice  of  expansion  and  testing  functions  al- 
lows efficient  calculation  of  matrix  elements  and  yields  accurate 
results  for  certain  widths  of  the  testing  function.  However,  for  other 
widths,  although  the  standard  criteria  for  the  selection  of  these 
(mictions  are  satisfied,  erroneous  results  are  obtained.  The  validity 
ol  the  moment  method  solution  can  be  checked  by  examining  the 
near-field. 
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New  Waveguide  Structures  for  Millimeter  Wave  and  Optical 
Integrated  Circuits 

ABSTRACT 

Some  new  dielectric  waveguide  structures  suitable  for  millimeter  wave 
and  optical  integrated  circuits  are  presented.  A method  of  analyzing  wave 
propagation  in  these  guides  is  developed  by  assuming  simple  field 
distribution  and  .approxiwal  i»g  the  various  regions  of  the  guides  in  terms  of 
effective  dielectric  constants.  The  mathematical  formulation  utilized 
results  in  simple  eigenvalue  equations  from  which  the  dispersion  character- 
istics of  Lhc  waveguides  arc  readily  obtained.  Experimental  results  arc 
described  and  the  agreement  between  theory  and  experiment  is  shown  to  be 
quite  good. 


Thl:;  work  was  supported  In  part  by  Joint  Service:;  Electronic?!  Program 
DAAliC- 0239  and  in  part  by  U.8.  Army  Research  Grant  DAWCO'i-'M-COI  1 3. 
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ANTENNAS  ON  COMPLEX  STRUCTURES 


Raj  Ml t era 

Electromagnetics  Laboratory 
Department  of  Electrical  Engineering 
University  of  Illinois 
Urbana,  Illinois  61801 


SUMMARY 

One  of  the  most  challenging  problems  facing  the  numerical  and  analyt- 
ical electromagnoticist  today  is  the  design  of  antennas  to  bo  mounted  on 
complex  structures,  such  as  automobiles,  ships  and  aircrafts.  The  problem 
is  complicated  for  at  least  two  reasons.  First,  the  modeling  of  complex 
structures  is  in  itself  a difficult  problem.  Second,  and  perhaps  the  more 
important  from  the  user's  point  of  view,  is  the  choice  of  the  type  of 
formulation  best  suited  for  the  particular  geometry  under  consideration. 
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Spectral  Theory  of  Diffraction* 

R.  Mittra,  Y Ruhmiit-Stimii.  and  W.  L.  Ko 
Electromagnetics  Laboratory 
Department  of  Electrical  Engineering 
University  of  Illinois 
Urbana.  III.  6IK0I,  USA 

Received  5 November  1975, Accepted  2 January  1976 

Abstract.  In  this  paper  we  present  an  overview  of  the  spectral  domain  approach  for  solving 
a variety  of  high  frequency  diffraction  problems.  We  demonstrate  via  a number  of  examples 
that  the  fields  derived  from  the  use  of  the  Spectral  Theory  of  Diffraction  (STD)  remain 
uniformly  valid  for  all  observation  angles,  including  the  shadow  boundary,  the  reflection 
boundary  and  the  caustic  direction  where  the  conventional  Geometrical  Theory  of  Diffrac- 
tion (GTD)  breaks  down. 

Furthermore,  we  show  how  the  accuracy  of  STD  can  be  improved  and  its  range  extended 
by  combining  it  with  integral  equation  methods  in  the  spectral  domain 

PACS  Codes:  84, 42,  10 
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Forward  Scattering  from  Square  Cylinders  in  the 
Resonance  Region  with  Application  to 
Aperture  Blockage 

w.  V.  T.  RUSCM,  SENIOR  mi.mm.r,  ii.rL,  JORGHN  APPIIL-I I ANSEN,  memoir,  n:rn,  CHARLES  A.  KLEIN, 

ANI)  RAJ  MITTRA,  FELLOW,  JLLC 


Abstract — The  relationship  between  the  induced  field  ratio  OH*)  of  a 
cylinder  and  aperture  blocking  of  a constant-phase  apn  t urc  by  c>  lindrical 
struts  is  discussed.  An  analytical  technique  is  presented  whereby  the  H R 
of  rectangular  c)lindcrs  can  be  calculated  using  the  mctliod-of- -moments 
with  internal  constraint  points.  An  experimental  technique  using  a 
forward-scatter  inp  range  is  used  to  measure  the  I Fit's  of  square  and 
circular  cylinders  in  an  anccltoic  chamber.  These  experimental  results  are 
compared  with  the  theory,  and  their  implications  on  aperture  blocking 
losses  and  borcsight  cross  polarization  are  discussed. 


Digital  Imaging  of  Gamma-Ray  Sources  with  Depth  Information 

KUSSIiLL  M SINGLLTON,  student  memiiek,  iiir.E,  PRHSTON  !..  RANSOM,  mi.muf.k,  ii:ee,  and 

RAJ  MITTKA,  i ellow,  ilee 


Abstract • A digital  method  for  imaging  gamma-ray  emitting  organs 
recorded  through  a Fresnel  rone  plate  (FZI>)  aperture  is  described. 
This  technique  can  be  used  in  an  almost  real-time  recording  and  imaging 
system  to  yield  quantitative  and  qualitative  cross  section  and  depth  in- 
formation of  the  source.  Furthermore,  this  method  avoids  some  of  the 
problems  introduced  hy  the  usual  method  of  optical  reconstruction. 
Experimental  results  of  the  recording  of  a multiplane  gamma-ray  source 
and  the  reconstruction  of  digitally  and  optically  formed  images  arc  also 
presented. 
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IMACE  RECONSTRUCTION  OF  MULTIPLANE  GAMMA-RAY  SOURCES 
WITH  A FRESNEL  ZONE  PLATE  APERTURE 


Raj  Mittra,  Russell  M.  Singleton,  Preston  L.  Ransom 
Electromagnetics  Labora tory 
Department  of  Electrical  Engineering 
University  of  Illinois 
Urbana,  Illinois  61801 

ABSTRACT 

A process  is  described  for  digitally  reconstructing  images  of  gamma- 
ray  emitting  organs  recorded  through  an  off-axis  Fresnel  zone  plate 
aperture.  This  technique  produces  image  cross  sections  for  given  depths 
which  contain  both  quantitative  and  qualitative  information  about  the 
source.  The  practicality  of  digital  decoding  is  demonstrated  experimentally 
by  comparing  the  digital  and  optical  reconstructions  of  images  of  a multi- 
plane gamma-ray  source.  Specifically,  the  multiplane  source  consists  of 
numerals  one,  two  and  three  located  in  three  distinct  planes  at  successive 
depths  from  the  detector.  The  digital  method  of  reconstruction  has  the 
following  advantages  over  a method  employing  optical  reconstruction.  A 
quasi-rcal  time  system  capable  of  recording  and  producing  images  can  be 
assembled  using  components  currently  available.  Quantitative  data  such 
as  the  source  size  and  depth  are  readily  determined.  Digital  image  process- 
ing can  be  easily  implemented  to  reduce  the  effect  of  artifacts  and  noise. 
Unlike  the  case  using  optical  reconstruction,  the  imago  depth  dimension  is 
not  distorted  In  comparison  with  the  source  depth.  Nonlinearity  and  film 
noise  associated  with  the  photographic  process  of  obtaining  a reduced  trans- 
parency are  avoided. 
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A TRANSFORM  APPROACH  TO  ELECTROMAGNET I C SCATTKR I NC  PROBLEMS 


Raj  Mittra 

El ec t romagnct  1 cs  l.abo rater y 
University  of  Illinois 
Urbana,  Illinois 


Generally  speaking,  there  are  two  distinctly  different  approaches 
available  for  solving  electromagnetic  scattering  problems.  Asymptotic 
methods,  for  large  wavenumber  k,  introduced  by  Keller  [1]  and  developed 
by  a great  many  researchers  [2,3],  are  usually  employed  for  scatterers 
whose  dimensions  are  large  in  terms  of  the  wavelength  of  the  incident 
^ield.  In  this  analysis,  scattering  is  considered  a local  phenomenon, 
identifiable  with  specific  parts  of  the  object,  e.g.,  specular  points, 
corners,  and  shadow  boundaries.  The  diffracted  far  fields  are  expressed 
in  terms  of  rays  for  which  the  diffraction  and  launching  coefficients  are 
determined  by  solving  canonical  problems  that  closely  resemble  the  surface 
configuration  in  the  vicinity  of  the  point  from  which  the  diffracted  field 
originates.  The  second  approach,  based  on  an  integral  equation  formulation 
for  the  unknown  surface  current  on  the  scatterer,  is  a global  formulation 
in  which  the  induced  current  distribution  on  the  entire  surface  of  the 
scatterer  must  be  solved  for  before  the  far  field  can  be  constructed.  The 
integral  equation  has  the  great  advantage  of  being  entirely  self-consistent 
as  the  boundary  condition  is  already  built-in  in  this  approach.  The  method 
of  solution  of  the  integral  equation  is  typically  based  on  the  moment  method 
which  transforms  the  original  equation  into  a matrix  form  that  is  numerically 
inverted  on  the  computer.  Though,  in  principle,  the  integral  equation  tech- 
nique is  valid  for  the  entire  frequency  range,  its  application  is  restricted 
to  frequencies  below  the  so-called  "resonance  region,"  i.e.,  when  the  dimen- 
sions of  the  scatterer  are  on  the  order  of  a wavelength  or  less.  This 
restriction  is  necessary  because  the  matrix  method  of  solution  becomes  very 
expensive  and  unwieldy  as  the  frequency  is  allowed  to  approach  the  resonance 
region.  On  the  other  hand,  high  frequency  asymptotic  methods  typically  be- 
come quite  inaccurate  when  applied  to  the  lower  frequency  range.  Furthermore, 
they  provide  no  convenient  accuracy  checks,  and  the  systematic  improvement  of 
these  solutions  may  also  be  difficult.  Thus,  there  exists  a great  need  for  a 
method  which  will  not  only  bridge  the  existing  gap  between  the  low  and  high 
frequency  regions,  but  will  also  provide  an  accuracy  check  and  a systematic 
improvement  procedure  for  the  high  asymptotic  solutions. 

The  purpose  of  this  paper  is  to  demonstrate  that  a recently  developed 
[3]  transform  method  of  solution  of  t fie  integral  equation  appears  to  meet  all 
of  the  desired  requirements  alluded  to  above.  The  application  of  the  Fourier 
transform  technique  to  the  solution  of  scattering  problems  is  by  no  means 
new.  The  problem  of  arbitrary  incident  waves  diffracting  from  a scatterer, 
c.g.,  a half-plane  lias  been  formulated  by  Cleinmnw  [A],  Felsen  and  Marcuvitz 
| 3]  and  Khestanov  [6],  to  name  just  a few.  The  Fourier  transform  technique 
itself  has  been  employed  by  Bojarski  [7]  for  solving  scattering  problems. 


QUASI-OPTTCAL  MILLIMETER  WAVEGUIDES  AND  COMPONENTS 


Raj  Mlctra 


ABSTRACT 

In  this  paper  we  review  a number  of  qunsi-optical  waveguide  struc- 
tures and  discuss  their  application  to  millimeter-wave  integrated  circuits. 
Some  novel  designs  for  dielectric  waveguides  and  passive  components  using 
these  waveguides  are  described.  A summary  of  the  results  of  theoretical 
and  experimental  studies  of  these  waveguides  and  circuits  is  presented. 
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A NFW  l.OOK  AT  TIIK  1 NTHCKAL  EQUATION  SOLUTION  OF  111011  FKFQIJKNCY  DIFFRACTION 
PROBLEMS 


Raj  Mittra 


ABSTRACT 

In  this  paper  we  report  a transform  method  for  eomhininp,  the  integral 
equation  and  high  frequency  asymptotic  techniques,  e.g.,  the  goomet r i ca  1 
theory  of  diffraction  or  GTD.  The  method  takes  advantage  of  the  fact  that 
the  Fourier  transform  of  the  unknown  surface  current  distribution  is  pro- 
portional to  the  scattered  far  field.  Two  methods  are  developed  for  sys- 
tematically improving  the  initial  form  of  the  high  frequency  asymptotic 
solution  by  manipulating  the  integral  equation  in  the  Fourier  transform 
domain.  Two  salient  features  of  the  transform  method  are  that  it  provides 
a convenient  validity  check  of  the  solution  and  that  it  yields  both  the 
induced  surface  current  density  as  well  as  the  far  field.  Several  illus- 
trative examples  that  demonstrate  the  usefulness  of  the  approach  for 
handling  a variety  of  electromagnetic  scattering  problems  in  the  resonance 
region  and  above,  and  comparison  with  other  methods  are  included  in  this 
p-.per. 
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Passive  Millimeter-Wave  IC  Components  Made  of 
Inverted  Strip  Dielectric  Waveguides 

**  R-  RUDOKAS  ano  T.  ITOH,  mnio*  mcmiuk,  ieie 

Ahitrcct — New  ilircction.il  couplers  :iml  riu"  resonators  for  niilliinctcr- 
wave  IC’s  were  fabricated  from  the  inverted  strip  (IS)  dielectric  wave- 
guide. Tlicy  were  tested  in  the  75  SO-GI 1/  ranyc,  and  the  agreement 
between  the  theoretical  and  experimental  results  was  found  to  be  good. 


Manuscript  received  March  17.  1076;  revised  July  51.  )07C.  This  work 
wav  vupi’tiiicd  in  I'.iil  I'V  the  Jiiuri  Services  I Icctroiiics  l*r.i|T.iin  I > A A 1 IU 7- 
iC,<.,,'!l,l,!v',n,‘l.',n,  f'"1  by  ,htf  U v A,n’)'  Research  Ullice  under  Grant 

tJ  A 111  HI-  .01  I .1, 

It  Kndol. ,iv  IV  vvnli  %tic  Coordinated  Science  t.ahorjiory  and  the  IVpart- 
"*•;»  •"  I lvvtrie.il  l.npneerini;.  Uimerviiyol  lllimuvji  Uilun.i  t liuinnairo 
Urhjii.i.  II.  IiIMII.  ' *• 

I.  Iioli  W'.iv  Will!  the  Coordinated  Svicncc  I ahoraiorv  and  ihc  Depart- 
ment  0|  | lev  I in  at  I iwmccruii:.  Uniwrilv  u|  llliiuiiv.il  lirli.iii.i-<  li  .mr.iirn 
I Ji liana,  II.  n|MI|.  lie  is  nmv  wnli  the  Manlurd  Kevealvh  liivlilule.  Menlo 
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ILIC  TKAK'.ACTIONS  on  MIC KOWAVE  HII'OKY  ANI)  IKIINKjULI,  VO|..  *07-24,  NO.  II,  NOVI  Mill  H I97f> 


Inverted  Strip  Dielectric  Waveguide  for 
Millimeter-Wave  Integrated  Circuits 

TATSUO  noil,  siniok  mcmulr,  inn 


Abstract — A new  type  of  dielectric  wavepuide,  which  has  a number  of 
advantages  out  other  previously  available  wavepuide  structures  for 
millimeter-wave  integrated  circuits,  is  described.  Dispersion  character- 
istics and  the  held  distributions  in  the  wavepuide  are  calculated  usinp  the 
concept  of  effective  dielectric  constant.  I ield  distributions  hasc  been 
measured  in  the  K0-(«ll/.  ranj'c  in  order  to  cluck  the  accuracy  of  the 
•nal)tical  results.  This  measurement  has  b en  done  usinp  a novel  experi- 
mental tcchni(|uc,  which  should  also  he  applicable  to  many  other  milli- 
meter-wave waveguides  and  components. 
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ABSTRACT 


Numerical  techniques  for  solutions  to  the  problems  of  electromagnetic 
radiation  and  scattering  are  considered  for  a finite,  hollow,  circular 
cylinder  of  radius  a.  The  singular-integral  equations  of  electromagnetic 
scattering  theory  arc  derived  along  with  their  extensions  to  thin 
surfaces  and  surfaces  with  edges.  In  addition,  constraints  are 
presented  which  are  necessary  for  a unique  solution  to  the  scattering 
problems  of  thin  structures.  The  equations  for  a finite  hollow  cylinder 
are  obtained  by  expanding  the  field  quantities  in  Fourier  series  about 
the  cylinder  axis  giving  rise  to  a separate  set  of  singular  integral 
equations  for  each  harmonic. 

The  method  of  moments  is  presented  as  the  basic  technique  of 
digitizing  the  integral  equations  for  numerical  solution.  It  is  found 
that  the  variational  interpretation  of  the  method  of  moments  can  be  used 
as  a guide  for  choosing  the  basis  and  testing  functions.  Of  particular 
interest  as  basis  functions  are  the  spline  functions  of  finite  support. 

The  spline  function  properties  of  smoothness  and  best  fit  are  also 
presented. 

It  is  shown  that  Hallen's  and  l’ockl ington' s formulations  for  the 
thin  wire  problem  are  equivalent  numerically  for  appropriate  testing 
functions.  It  is  also  shown  that  the  I’ocklington  form  is  more  desirable 
when  smooth  basis  functions  are  used  in  conjunction  with  pulse  or  delta 
testing  functions.  In  this  context,  the  second-order  sinusoidal  spline 
is  found  to  be  an  excellent  current  representation  for  both  the  scattering 
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and  pulse-feed  problems.  Slope  dlscontluull les  are  easily  included  to 
approximate  a delta-feed  problem.  Approximate  operators  are  also 
considered  with  emphasis  on  the  equivalence  of  the  finite  difference 
approximation  and  piecewise  sinusoidal  basis  functions. 

The  problems  of  coupling  in  the  higher  harmonics  are  investigated 
for  the  first-harmonic  problem.  It  is  shown  that,  minimum  coupling  of 
the  equations  is  desirable  in  addition  to  the  dominance  of  each  equation 
over  the  entire  structure  by  its  respective  current  component.  These 
features  are  obtained  using  a new  set  of  equations  obtained  from  the 
combination  of  the  equations  for  the  tangential  electric  field  and  the 
normal  magnetic  field.  These  equations  are  related  to  the  normal 
derivative  of  the  tangential  magnetic  field  equations  which  are 
well-behaved  for  thick  structures.  The  solutions  of  these  equations 
are  in  excellent  agreement  with  the  results  of  other  workers.  Various 
operator  approximations  are  also  considered. 
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ABSTRACT 

An  efficient  method  ha9  been  developed  for  performing  the 
Karhunen-Loeve  transform  which  is  optimum  for  reducing  the  number  of 
samples  to  be  sent  via  a given  channel  for  high  fidelity  picture 
transmission.  The  method  is  based  on  the  use  of  the  Fast  Fourier 
Transform  algorithm  and  linear  interpolation  procedure.  Some  test 
computations  illustrating  the  use  of  the  method  are  included  in  the 
paper. 
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EXCITATION  OF  AN  OPTICAL  FIBER  BY 
A CALL; SIAN  BEAM 


* + * 

M.  Mostafavi  , T.  Itoh  and  R.  Mittra 

*Department  of  Electrical  Engineering 
University  of  Illinois,  Urbana,  Illinois  61801 
Department  of  Electrical  Engineering  and  Coordinated  Science  Laboratory 
University  of  Illinois,  Urbana,  Illinois  61801 


ABSTRACT 

In  this  paper  we  present  t lie  results  of  an  in-depth  analysis  of  the 
problem  of  Gaussian  beam  excitation  of  optical  fibers  at  normal  incidence. 

The  first  step  involves  t lie  derivation  of  the  eigenvalues  and  eigenfunc- 
tions of  the  guided  modes  as  well  as  the  radiation  field  of  the  fiber.  Then, 
the  eigenvalue  equation  is  solved  for  the  frequencies  of  interest  using 
numerical  techniques.  The  second  step  entails  the  application  of  the  con- 
tinuity conditions  at  the  interface  located  at  the  fiber  end  using  Gaussian 
beams  with  different  spot  sizes  as  our  excitation  source.  Next,  the  simul- 
taneous set  of  equations  obtained  is  numerically  solved  for  the  reflected, 
guided  and  radiated  modal  coefficients  and  their  corresponding  powers. 


This  work  was  supported  by  the  U.  S.  Army  Research  Office  tinder  Grant  No. 
DAIIC0W/.-C-0113. 
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ABSTRACT 


The  method  of  moments  is  a simple  technique  for  converting  integral 
equations  occurring  in  electromagnetics  problems  into  matrix  equations 
which  can  be  easily  solved  on  digital  computers.  Although  the  method  of 
moments  can  be  readily  applied  to  a wide  class  of  problems,  it  should 
not  be  assumed  that  the  computed  results  automatically  converge  to  the 
correct  solution.  One  of  the  difficulties  that  can  occur  is  an  instability 
of  the  method  to  small,  unavoidable  errors.  To  measure  this  possible 
ill-conditioning,  stability  parameters  are  discussed.  The  primary 
parameter  used  will  be  the  matrix  condition  number,  a scalar  function  of 
the  matrix,  which  is  useful  in  estimating  errors  in  the  calculated 
solution  due  to  errors  in  either  the  matrix  or  the  excitation  vector. 

The  a posteriori  condition  number  provides  an  improved  estimate  by 
considering  the  particular  excitation  vector.  The  pivot  ratio  has 
definite  computational  advantages.  Several  methods  for  improving  the 
stability  of  the  matrix  equation  are  compared. 

The  discussion  and  theory  of  stability  and  convergence  are  applied 
to  six  areas  of  applications  — waveguide  discontinuities,  wavefront 
reconstruction,  electrostatics,  remote  sensing,  scattering  by  rectangular 
conducting  cylinders,  and  thin-wire  antennas  and  scatterers.  The  chapters 
on  these  applications  not  only  demonstrate  the  theory  but  also  discuss 
other  investigations  of  the  method-of-moments  technique.  The  best  method 
for  calculating  Induced  Field  Ratios,  a useful  parameter  in  describing 
scattering  by  rectangular  cylinders,  is  presented.  A new  application  of 
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the  method  of  moments  to  thin-wire  modeling  of  antennas,  which  is  using 
piecewise  sinusoids  and  rectangular  pulses  as  modeling  functions,  is 
compared  with  the  reaction-matching  technique.  The  surprising  equivalence 
between  finite  differencing  and  the  use  of  sinusoids  with  pulses  is 
discussed.  Also  studied  are  antennas  with  junctions  and  straight 
antennas  modeled  with  unequal  length  segments. 
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ABSTRACT 


The  finite  length  digital  deconvolution  problem  Is  formulated  and 
discussed  in  terms  of  modern  optimization  theory.  The  ill-conditioned 
nature  of  deconvolution  is  identified  and  classic  deconvolution  operators 
are  examined  in  terms  of  their  respective  effectiveness  and  ease  of  imple- 
mentation. 

Continuous  deconvolution  filters  that  are  optimal  in  terms  of  min- 
imum variance  and  solution  smoothness  are  derived.  These  filters  contain 
digital  counterparts  which  can  be  obtained  by  means  of  conventional  fre- 
quency sampling  techniques. 

Quantization  and  frequency-aliasing  are  identified  as  characteristic 
noise  sources  of  digital  deconvolution  and  analyzed  in  terms  of  their  re- 
spective corrupting  effects  upon  accurate  optical  signal  estimation. 

Three  approaches  for  eliminating  frequency-aliasing  error  arc  discussed. 

A multiple  intensity  measurement  technique  for  determining  the  magnitude 
and  phase  of  complex  optical  wavefronts  by  means  of  a coherent  optics-mi- 
croprocessor hybrid  system  is  developed. 

A picture-shift  sampling  scheme  is  proposed  for  effectively  increasing 
picture  sampling  resolution  without  requiring  additional  computer  core.  A 
method  is  developed,  based  upon  this  scheme  and  incoherent  light  illumi- 
nation, for  performing  continuous  low-pass  filtering.  A noise  analysis 
of  the  sampling  technique  is  later  performed. 

linear  system  testing  is  discussed  in  terms  of  its  role  with  respect 
to  digital  deconvolution  processing.  Test  accuracy  requirements  are  for- 
mulated to  guarantee  the  possibility  of  success  I u.l  deconvolution  algorithm 
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implementation.  An  optimal  test  signal  and  testing  strategy  for  linear 
system  measurement  are  analytically  derived  by  maximizing  the  probability 
of  acceptable  deconvolution  estimation. 
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ABSTRACT 


This  paper  introduces  a new  approach  for  combining  the  integral 
equation  and  high  frequency  asymptotic  techniques,  e.g.,  the  geometrical 
theory  of  diffraction.  The  method  takes  advantage  of  the  fact  that  the 
Fourier  transform  of  the  unknown  surface  current  distribution  is  pro- 
portional to  the  scattered  far  field.  A number  of  asymptotic  methods 
are  currently  available  that  provide  good  approximation  to  this  far 
field  in  a convenient  analytic  form  which  is  useful  for  deriving  an 
initial  estimate  of  the  Fourier  transform  of  the  current  distribution. 

An  iterative  scheme  is  developed  for  systematically  improving  the 
initial  form  of  the  high  frequency  asymptotic  solution  by  manipulating 
the  integral  equation  in  the  Fourier  transform  domain.  j 

A synthetic-aperture-distribution  scheme  is  also  developed  in  which 
the  approximate  scattered  far-field  pattern  obtained  by  asymptotic  tech- 
niques is  improved  by  systematically  correcting  the  scattered  field 
distribution  on  an  aperture  erected  in  juxtaposition  with  the  obstacle. 
The  introduction  of  such  a planar  aperture  not  only  provides  an  additional 
degree  of  freedom  in  performing  improving  operations,  but  also  renders 
the  scheme  to  handle  n-diincnsional  geometries  by  (n  - l)-dimensional 
fast  Fourier  transform  (FFT),  where  n = 2,3,  and  circumvents  the  un- 
wieldy three-dimensional  FFT,  making  it  a conceptually  simple  and  com- 
putationally efficient  method. 
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ABSTRACT 


The  knowledge  of  high- f rcqucncy  di f f met  ion  of  an  arbitrary  wave  inci- 
dent on  an  edge  is  important  in  many  applications,  such  as  antennas  mounted 
on  aircraft  and  reflector  antennas  illuminated  by  complex  feeds.  In  this 
paper  the  problem  of  a half-plane  illuminated  by  a nonplanar  wave  is  investi- 
gated using  the  concept  of  the  plane  wave  spectral  representation.  Tor  large 

wave  number  k,  a new  higher-order  asymptotic  solution  for  the  total  field 

-5/2 

up  to  and  including  terms  of  order  k relative  to  the  incident  field  is 

derived.  The  behavior  of  the  solution  for  the  observation  points  which 
coincide  with  shadow  boundary  directions  of  a multipole  line  source  is 
discussed  in  detail.  Furthermore,  numerical  solution  of  the  field  integral 
representation  is  constructed  for  the  observation  angles  in  the  transition 
regions.  The  results  are  compared  with  those  of  the  Geometrical  Theory  of 
Diffraction  (GTD) , the  Uniform  Asymptotic  Theory  (UAT) , the  Uniform  Theory 
of  Diffraction  (UTD)  and  the  Modified  Slope  Diffraction  (MSD) . 
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A SPECTRAI,  DOMAIN  INTERPRETATION  OF  MICH 
FREQUENCY  DIFFRACTION  PHENOMENA 


Y.  Rahmat-Samii  and  R.  Mittra 
DSpattment  of  Electrical  Engineering 
iTfrivcrsity  of  Illinois 
Urbana /'Illinois  61801 


ABSTRACT 

A spectra]  domain  interpretation  of  high  frequency  diffraction  phenomena 
is  discussed  introducing  the  concept  of  a spectral  diffraction  coefficient 
which  resembles  Keller’s  coefficient.  The  solution  of  the  two-dimensional 
problems  of  diffraction  of  an  arbitrary  field  (with  no  caustics)  by  a half- 
plane is  investigated  and  results  are  given  for  any  observation  angle  including 
in  particular,  the  determination  of  the  field  at  the  shadow  boundaries.  The 
high  frequency  scalar  diffraction  by  apertures  and  semi- inf  ini te  cylinders  is 
formulated  in  a systematic  manner  and  the  formulation,  which  is  valid  for  any 
observation  angle,  is  compared  with  that  of  Ufimtsev's.  Results  are  also 
given  for  the  diffracted  field  at  the  caustics. 


ON  THE  INVESTIGATION  OK  DIFFRACTED  FIELDS  AT 
THE  SHADOW  UOUNDAK  1 ES  OK  STAGGERED  J’AKALLEL 
l'LATLS— A Sl’ECTKAL  DOMAIN  A l’ 1’ KOACII 


Y.  Rahmal-Samii  and  K.  Mittra 
Electromagnet  ics  Laboratory 
Department  of  Electrical  Engineering 
University  of  Illinois  at  Urbana-Cliampu  i gn 
Urbana,  Illinois  61801 


ABSTRACT 

The  diffraction  of  plane  waves  by  two  staggered  parallel  plates 
is  investigated  using  the  spectral  diffraction  coefficient  in  the  Fourier 
transform  domain.  It  is  shown  that  when  the  interaction  between  the 
plates  is. negligible , the  total  field  can  be  represented  in  terms  of  a 
double  Fourier  integral,  a result  that  is  identical  to  the  one  derived 
by  Jones  via  the  Woiner-Hopf  technique.  The  double  integral  is  asymp- 
totically evaluated  via  the  saddle-point  integration  technique  and  com- 
pact expressions  for  the  field  are  derived,  up  to  the  order  k \ .for 
observation  points  located  at  the  shadow  boundaries  and  also  away  from 
the  transition  regions.  These  expressions,  being  exact,  serve  a very 
useful  purpose,  viz . , they  are  used  to  test  different  uniform  and 
nonuniform  asymptotic  theories  at  the  shadow  boundaries.  The  results 
of  these  tests  are  included  in  the  paper. 
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A SPECTRAL  DOMAIN  METHOD 
FOR  REMOTELY  PROBING 
STRATIFIED  MEDIA 


ABSTRACT 

The  problem  of  remotely  probing  a strutified,  lossless,  dielectric  medium 
is  formulated  using  the  spectral  domain  method  of  probing.  The  response  of 
the  medium  to  a spectrum  of  plane  waves  incident  at  various  angles  is  used  to 
invert  the  unknown  profile.  For  TE  polarization,  the  electric  field  satisfies 
a Helmholtz  equation.  The  inverse  problem  is  solved  by  means  of  a new  repre- 
sentation for  the  vrave  function.  The  principal  step  in  this  inversion  is  solv- 
ing a second  kind  Fredholm  equation  which  is  very  amenable  to  numerical  computa- 
tions. Several  examples  are  presented  including  some  which  indicate  that  the 
method  can  be  used  with  experimentally  obtained  data.  When  the  fields  exhibit 
a surface  wave  behavior,  a unique  inversion  can  be  obtained  only  if  information 
about  the  magnetic  field  is  also  available.  In  this  case,  the  inversion  is 
accomplished  by  a tvo-step  procedure  which  employs  a formula  of  Jost  and  Kohn. 
Some  examples  are  presented  and  an  approach  which  greatly  shortens  the  compu- 
tations without  greatly  deteriorating  the  results  is  discussed. 


TRANSMISSION  OF  AN  111  WAVE  THROUGH  THE  APERTURE 
OF  A CYLINDRICAL  CAVITY* 


S.  Safavi-Naini , S.  W.  Lee  and  R.  Mittra 
Electromagnetics  Laboratory 
Department  of  Electrical  Engineering 
University  of  Illinois 
Urbatia,  Illinois  61801 


ABSTRACT 

The  penetration  of  an  electromagnetic  wave  through  an  aperture 
in  a cylindrical  structure  is  investigated.  By  using  the  moment  method, 
the  electrical  and  magnetic  fields  in  a cavity  behind  the  aperture  are 
determined  as  a function  of  frequ  ncy  of  the  incident  field.  The  accuracy 
of  the  numerical  solution  is  established  through  tests  of  satisfaction 
of  the  boundary  condition,  edge  condition,  and  convergence  with  respect 
to  the  number  of  modes  in  the  aperture.  Depending  on  the  frequency  and 
location  inside  the  cavity,  the  stored  electromagnetic  energy  density 
varies  very  rapidly.  Its  peak  value  can  be  two  orders  of  magnitude 
greater  than  the  incident  energy  density.  The  frequencies  where  the 
peaks  occur  can  be  identified  approximately  as  the  resonance  frequencies 
of  the  cavity  when  the  aperture  is  closed. 
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AN  APPROACH  TO  HIGH-FREQUENCY  SCATTERING  FROM  SMOOTH  CONVEX  SURFACES 


By 

R.  Mittra  and  W.  L.  Ko,1[ 

F.lectromagnct  ics  Laboratory 
Department  of  Electrical  Engineering 
University  of  Illinois  at  Urbnna-Champaign 
Urbana,  Illinois  61801 

ABSTRACT 

In  this  paper  we  introduce  a method  based  on  the  Galerkin  procedure 
in  the  transform  domain  that  not  only  allows  accurate  solutions  for 
both  the  far-field  and  the  induced  surface  current  be  derived  but  also 
provides  a convenient  accuracy  check  for  the  solutions  thus  obtained. 

The  total  solution  is  constructed  using  only  the  geometrical  optics 
method  and  a newly  introduced  synthetic  aperture  technique  well  suited 
for  deriving  efficient  transform  solutions  for  n-dimensional , smooth 
curved  objects  using  (n  - 1)  - dimensional  transforms.  The  need  for  the 
conventional  creeping-wave  contribution  is  obviated  in  this  method. 


This  work  was  supported  by  the  U.S.  Army  Research  Office  under  Grant  No. 
DAHC04-74-G-01 1 3 and  the  Office  of  Naval  Research  under  Grant  No. 

N000 1 4-7  5-C-0293 . 
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A NEW  APPROACH  TO  THE  111  IN  SCATTERED 
PROBLEM  USING  THE  HYBRID  EQUATIONS 

William  A.  Davis,  Member,  IEEE,  and  Raj  Mittra,  Fellow,  IEEE 

ABSTRACT 

A new  set  of  integral  equations  for  electromagnetic  scattering 
problems,  the  "hybrid"  equations,  are  presented.  The  advantages  of 
these  equations  for  thin  perfect  conductors  are  discussed  in  com- 
parison to  the  magnetic  and  electric  field  integral  equations. 
Specific  comparisons  are  made  with  the  solution  of  the  electric 
field  integral  equation  for  a finite  hollow  cylinder.  It  is  demon- 
strated that  the  primary  advantage  of  these  equations  is  obtained 
by  minimizing  the  coupling  between  component  equations  for  the  two 
surface  currents. 
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A New  Approach  Based  on  a Combination  of  Integral 
Equation  and  Asymptotic  Techniques  for  Solving 
Electromagnetic  Scattering  Problems 
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Abilracl -Wf  introduce  a new  approach  for  combining  die  integral 
equation  and  high  frequency  asymptotic  technique',  e.g..  the  geometri- 
cal theory  of  (Infraction.  The  method  lakes  advantage  of  the  fact  that 
the  I onrier  transform  of  thb  unknown  surface  current  dislrihiilion  is 
proportional  to  flic  scattered  f tie  - field . A number  of  asymptotic  meth- 
ods ate  currently  available  that  provide  pood  approximation  to  Ibis  fas- 
ficld  in  a convenient  analytic  form  vsUu  li  is  useful  for  deriving  an  initial 
estimate  of  the  I onrier  transform  of  lire-  current  distribution.  An  itera- 
tive scheme  is  developed  for  systematically,  improving  the  initial  form 
of  the  high  frequency  asymptotic  solution  hybiji  aiiipiil.iling  the  integral 
equation  in  the  Fourier  transform  domain.  A .salient  feature  of  the 
method  is  that  it  provides  a convenient  validity  check  of  the  solution 
for  the  surface  current  distribution  by  verifying  that  ijie  scattered  field 
il  ladialcs  indeed  satisfies  the  boundary  conditions  at  tile  suif.icc  of  the 
scitlercr.  Another  important  feature  of  the  rnc-l hod  is  that  it  yields 
both  the  induced  surface  current  density  anil  tile  fat  field,  diffraction 
by  > strip  (two-dimensional  problem)  and  diffraction  by  a"Hjin  plate 
(threcdimension.il  problem)  are  presented  as  illuslralise  examples  that 
demonstrate  the  usefulness  of  the  approach  for  handling  a variety  of 
electromagnetic  scattering  problems  in  the  resonance  region  ami  abugc. 
Some  concluding  remarks  and  comparison  with  oilier  methods  arc  als<A 
included. 
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